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Abstract 

Climate forcing by anthropogenic aerosols is a significant but highly uncertain factor in global 
climate change. Only satellites can offer the global coverage essential to reducing this 
uncertainty; however, satellite measurements must be coupled with correlative, in situ 
measurements both to constrain the aerosol optical properties required in satellite retrieval 
algorithms and to provide chemical identification of aerosol sources. This grant funded the first 
two years of a three-year project which seeks to develop methodologies for combining 
spacebome lidar with in-situ aerosol data sets to improve estimates of direct aerosol climate 
forcing. Progress under this two-year grant consisted in the development and deployment of a 
new in-situ capability for measuring aerosol 180° backscatter and the extinction-to-backscatter 
ratio. This new measurement capacity allows definitive lidar/in-situ comparisons and improves 
our ability to interpret lidar data in terms of climatically relevant quantities such as the extinction 
coefficient and optical depth. Measurements were made along the coast of Washington State, in 
Central Illinois, over the Indian Ocean, and in the Central Pacific. Thus, this research, combined 
with previous measurements by others, is rapidly building toward a global data set of extinction- 
to-backscatter ratio for key aerosol types. Such information will be critical to interpreting lidar 
data from the upcoming PICASSO-CENA 1 , or P-C, satellite mission. Another aspect of this 
project is to investigate innovative ways to couple the lidar-satellite signal with targeted in-situ 
measurements toward a direct determination of aerosol forcing. This aspect is progressing in 


1 Pathfinder Instruments for Cloud and Aerosol Spacebome Observations - Climatologie Etendue des Nuages et des 

Aerosols 
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collaboration with NASA Langley's P-C lidar simulator and radiative transfer modeling by the 
University of Lille, France. 


Background and Rationale 

Progress on the problem of climate forcing by anthropogenic aerosols will require global- 
scale observations, which, in turn, can only be provided by satellite-borne instruments. With this 
in mind, a joint U.S./French research satellite project, PICASSO-CENA, or P-C, is being 
developed for a 3-year mission beginning in 2003. Following the recommendation of the U. S. 
National Research Council (NRC, 1996), this project will feature Light Detection and Ranging 
(lidar) instrumentation for obtaining quantitative, vertically resolved aerosol information both 
day and night, over both continents and oceans. The capabilities of spacebome lidar have been 
studied previously via deployment on the U. S. Space Shuttle (Winker et ah, 1996). 

Lidar has become a central technology in current strategies for tropospheric aerosol 
research because of its demonstrated ability to map aerosol variations throughout the atmospheric 
column. Its use is complicated, however, by the fact that the lidar signal contains a convolution 
of two basic optical properties of the aerosol particles: the 180° backscatter coefficient, (m 
sr ') and the extinction coefficient, G ep (m 1 ). A quantitative retrieval of either property requires 
knowledge their relationship along the laser path. A central goal of this project is to review and 
extend current knowledge of this relationship, which we express in the usual form of an 
extinction-to-backscatter ratio, S = G ep / (3 p , with units of sr. (To make this definition 
unambiguous, the value of S in the Rayleigh limit is 8n/3.) 

Like other aerosol optical properties, S can be calculated from Mie Theory. Mie 
calculations are essential for instrument validation experiments (e.g. Fig. 1) and provide 
guidance to field measurements by suggesting the range of values to be expected and what the 
controlling factors ought to be. A study by Ackermann (1998) illustrates the latter point. He 
used assumptions for aerosol size distribution and refractive index to show that S should vary 
from about 15 to 75 sr for tropospheric aerosols, with lower values associated with weakly 
absorbing coarse-mode particles (i.e. seasalt and mineral dust) and higher values associated with 
small and/or highly absorbing accumulation-mode particles. These Mie calculations were 
combined with models of aerosol hydration to develop functional relationships between S and 
ambient relative humidity for several aerosol types. 

While valuable, such calculations are inadequate in themselves for two reasons. First, 
Mie calculations are limited to simple, frequently unrealistic particle morphologies 
(homogeneous spheres, concentric spheres, and a few variants). Second, this "theoretical" 
approach to determining S can never be more accurate than the empirical knowledge of particle 
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size and refractive index (the latter from knowledge of particle chemistry) upon which the 
calculation is based. For these reasons, direct empirical determinations of S are required. 

Several methods have been developed to make this measurement; however, as shown in 
Table 1, below, they have been deployed to date in only a handful of investigations of 
tropospheric aerosols. Most regions of the Earth and many major aerosol types have no 
published data at all and a statistically significant database exists for only one location - Tucson, 
Arizona. This project seeks to expand that data base in preparation for the upcoming PICASSO- 
CENA satellite mission. The last four entries in Table 1 show the progress that has been made so 
far. 

Progress 

This project was originally conceived in terms of three tasks. Progress is therefore grouped 
under those original headings, although the direction of research has clearly changed in light of 
the initial findings. 

Task 1 Determining lidar ratios from existing in situ optical measurements 

This task was superceded by the progress under Task 2. Results from direct, calibrated 
measurements showed that existing, nephelometric data sets combined with Mie Theory do not 
offer an accurate means of predicting the lidar ratio. Therefore, our focus has shifted to 
expanding the data base of direct measurements. 

Task 2 Design of a new instrument for directly measuring 180-backscatter 

Under this task we designed, built, calibrated, and deployed a new, nephelometer-type 
instrument for performing direct, calibrated measurements of 180° backscatter and the 
extinction-to-backscatter ratio, S. Results have been reported in one published article (Doherty 
et al., 1999) and one article submitted for publication (Anderson et al., 2000). A summary of 
those, plus more recent, results in included in this section. 

Figure 1 shows the calibration of the device with laboratory particles of known size and 
refractive index (latex spheres). The continuous lines represent Mie Theory prediction. 

Measured points are all seen to fall within experimental uncertainties, indicated by horizontal and 
vertical error bars. Based on this well-calibrated methodology, measurements were made at four 
sites in different regions of the globe. Results are summarized in Table 1, which serves to 
document as well the paucity of previous measurements of S.. , 
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Task 3. Analyze 180-backscatter optical depth as an index of aerosol climate 
forcing 

This task has evolved into a collaborative project involving the lidar simulator model 
developed at NASA Langley as well as radiative transfer modeling at the University of Lille, 
France. Progress to date has involved developing a set of atmospheric aerosol "scenes" to be 
jointly studied by the lidar simulator and the radiative transfer models, analyzing preliminary 
results from simple tests, and designing future simulation experiments. Considerable progress on 
this task is expected in year three of the project, which is funded under a separate grant and will 
thus be the subject of a future report. 

A critical underpinning to these experiments is accurate knowledge of P-C noise 
uncertainty and how this effects the spatial resolution of aerosol backscatter data. The remainder 
of this section details progress made in that area. 

The effective spatial resolution of P-C lidar for aerosol studies is limited by instrumental 
noise, which has now been estimated via the lidar simulator. Results for the 532 nm from initial 
runs using two simple aerosol models imply that the noise is well behaved for clear-sky 
conditions and at modest levels of aerosol extinction. The equations below describe a simple 
model which appears to fit the simulator data quite well. 

Backscatter by aerosol particles, P a , is derived as total measured backscatter, P, minus the 
calculated molecular backscatter, P m , 

Pa=P"P m 0) 

Therefore, the uncertainty in P a , which I will denote 8P a , is, 

8P a = A /8p 2 +5p m 2 (2) 

Assuming good knowledge of molecular density, the uncertainty in P m is just due to the 
calibration uncertainty, 5C, 

O) 

We set 5C/C to 10% for this calculation. We use the simulator results for the elevated layers of 
"dust" and "dirty pollution" to estimate 5p, using a criterion of one standard deviation (70% 
confidence) of the simulated power for various averaging scales with respect to the grand average 
over the entire aerosol layer (3-5 km) and over all simulated "shots". To a very good 
approximation, the standard deviation of the simulated power scales with the square root of the 
signal, such that we can use the noise uncertainty for molecular backscatter as a reference point, 

(4> 

In addition, simulated noise scales inversely with the square root of n, the number of "bins" over 
which the signal is averaged. Here we use the word "bin" to refer to the highest spatial resolution 
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of a single shot at 30 m vertical resolution. The results in the lower troposphere are very nearly 
fit by. 


SPm _ 1 

Pm Vn 


(5) 


Combining (2) - (5) and dividing through by P a , 


5Pa _ Pm I 

|iJ_ + 

fscf 

Pa Pa 1 

n Pm 

[c J 


Figure 2 shows contours of 5(3/P a plotted as a function of n (on the x-axis) and the magnitude of 
P, (on the y-axis). The y-axis is actually P a scaled by molecular backscatter, since this is a more 
meaningful quantity to most people. Given knowledge of the extinction-to-backscatter ratio, S a , 
the magnitude of P a can be converted to an aerosol extinction value, a,. 


a a (km _1 ) = 0.00163^-S a 
Pm 


( 6 ) 


(Here, the aerosol extinction is implicitly referenced to standard temperature and pressure, since 
0.00163 km' 1 sr' 1 is the 532 nm molecular backscatter at STP.) The right-hand y-axis shows the 
aerosol extinction level for S,=60 sr' 1 , a typical value for pollution aerosol. For this value of S,, 
P/Pn=l corresponds to an extinction level of 0.1 km' 1 , a very high level of extinction. The green 
box on Figure 2 indicates the region where reasonable accuracy (5-20%) is achieved for 
reasonable aerosol extinction levels (0.01 to 0.1 km' 1 ). According to these {very preliminary) 
calculations, averaging over 100 to 1000 "bins" will be necessary for accurate measurements of 
pollution aerosol backscatter outside of intense plumes. To put this in context, 200-bin 
averaging would correspond to averaging over 150 m in the vertical (5-vertical bits) and 10 km 
in the horizontal (40-shots). 

Based on this knowledge of the effective spatial resolution of the P-C aerosol backscatter 
data, we are in the process of designing a correlative in-situ measurement strategy for the P-C 
mission. We expect to submit this strategy for publication - and thus for consideration by the 
worldwide aerosol measurement community - within the next few months. 
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Tables 


Table 1: Measurements of S for tropospheric aerosols 


Method 

Location 

N <‘> 

S(sr) 

slant-path lidar (b) 

mixed layer, west U.S. 
(Tucson, AZ) 

81 

8-75 

horizontal lidar (c) 

Netherlands 

10 

10-50 

horizontal lidar (d) 

marine surface layer, 
Australian coast 

10 

40-80 

multi-wavelength 
lidar (e) 

marine boundary layer 
(tropical Atlantic) 

l 

<30 


mixed layer over 
rainforest (S. Amer.) 

1 

43-60 


Saharan dust aloft 
(tropical Atlantic) 

2 

15-62 

backscatter-sonde (f) 

rural, arid SW U.S. 

12 

42 


rural western U.S. 
lower troposphere 

4 

15-30 


rural western U.S. 
upper troposphere 

4 

15-60 

space-bome lidar (g) 

smoke layers in SH 
upper troposphere 

7 

50-90 

Raman lidar (h) 

polluted lower tropo. 
over Northern Germany 

2 

20-40 

Raman lidar (i) 

polluted boundary layer 
over Leipzig, Germany 

1 

55-95 

Raman lidar (j) 

rural, lower tropo. over 
Oklahoma 

25 

20-80 

this work (k) 

Northwest Coastal U.S. 

3 

20-70 

this work (1) 

Indian Ocean 

50 

40-90 

this work (m) 

polluted, Central U.S. 

70 

30-71 

this work (n) 

Central Pacific 

50 

28-32 


a Approximate number of independent samples. 

b Determines effective S over mixed layer at 694 nm. (Spinhirne et al., 1980; Reagan et al., 1984; 1988) 
c 1064 nm laser. No information on location of instrument or type of aerosol investigated, (de Leeuw et al., 
1986) 

d 532 nm laser beam aimed 2 m over ocean surface from coastal site. Effect of waves not assessed. (Young et al., 


1993). 

e S at 600 nm is constrained by Mie calculations based on the measured wavelength variation of P p . (Sasano and 
Browell, 1989) 

f S reported at 690 nm. Method requires significant wavelength and angular adjustments, based on assumed size 
distributions and Mie calculations. (Rosen et al., 1997a, b) 
g S at 532 nm is constrained to values yielding physically plausible lidar retrievals under the assumption of 
constant S throughout upper troposphere. (Kent et al., 1998) 
h S at 308 nm determined from independent extinction measurements using nitrogen-Raman; one profile. 
(Ansmann et al., 1992) 

i S at 532 nm using nitrogen-Raman as above. (Muller et al., 1998) 

j S at 351 nm using nitrogen-Raman as above; profiles in lower troposphere over 8 nights (Ferrare et al., 1998) 
k S at 532 nm using the nephelometric technique discussed herein. S was ca. 20 during marine flow and 60-70 
during continental flow. (Doherty et al., 1999) 

1 As in k except for airborne measurements during the Indian Ocean Experiment, 
m As in k except for surface measurements in Central Illinois, 
n As in k except for surface measurements on the windward coast of Oahu. 








